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We report on the measurement of the heat capacity for an optically-trapped, 
strongly-interacting Fermi gas of atoms. In the experiments, a precise input of 
energy to the gas is followed by single-parameter thermometry. The thermom- 
etry determines a temperature parameter f from the best fit of a Thomas- 
Fermi distribution with a fixed Fermi radius to the spatial density of the cloud. 
At T — 0.33, we observe a transition between two patterns of behavior: For 
f = .33 — 2.15, we find that the heat capacity closely corresponds to that of a 
trapped normal Fermi gas of atoms with increased mass. At low temperatures 
T = 0.04 — 0.33, the heat capacity clearly deviates from normal Fermi gas 
behavior. 

Strongly-interacting, degenerate atomic Fermi gases (1) provide a paradigm for strong in- 
teractions in nature ©. Measurements of the interaction energy (0|5]|?]|5l) test predictions of 
universal interactions in nuclear matter (0001, as well as effective field theories of strong in- 
teractions ©. The anisotropic expansion observed for strongly-interacting Fermi gases © is 



analogous to the "elliptic flow" of a quark-gluon plasma (10). High temperature superfluidity 
has been predicted (|i 1 1 li 21 1731 17?l 1751 [7(511771) in strongly-interacting Fermi gases, which can be 
used to test theories of high temperature superconductivity (ITffl) . Microscopic evidence for high 
temperature superfluidity has been obtained in the condensation of preformed pairs U9\\20b and 
in radio frequency measurements of the pairing gap j21\\22h . Macroscopic evidence arises in 
anisotropic expansion ( 1 ) and in collective excitations d23H24ll231) . 

In superconductivity and superfluidity, measurements of the heat capacity (26) have played 
an exceptionally important role in determining phase transitions d27l) and in revealing the nature 
of the many-body quantum state of the system. We report on the measurement of the heat 
capacity for a strongly-interacting Fermi gas of 6 Li atoms, confined in an optical trap. Our 
experiments examine the fundamental thermodynamics of the gas. In the following, we first 
describe how the gas is prepared and our method for adding a precisely known energy to the 
gas. Then we discuss our technique of thermometry, which provides a monotonic temperature 
scale and a well-defined method for comparing experiment with predictions. 

We prepare a degenerate 50-50 mixture of the two lowest spin states of 6 Li atoms by forced 
evaporation in an ultrastable CO2 laser trap (l2fft as described previously 0. At a bias magnetic 
field of 840 G, just above the Feshbach resonance, the trap depth is lowered by a factor of ~ 580 
in a few seconds (171123b and then recompressed to 4.6% of the full trap depth in 1.0 s and held 
for 0.5 s to assure equilibrium. After a controlled amount of energy is added to the gas, as 
described below, the gas is allowed to thermalize for 0.1 s. Finally, the gas is released from the 
trap and imaged at 840 G to determine the number of atoms and the temperature parameter T. 
The column density is obtained by absorption imaging of the expanded cloud after 1 ms time 
of flight, using a two-level state-selective cycling transition (171 1231). In the measurements, we 
take optical saturation into account exactly and arrange to have very small optical pumping out 
of the two-level system. For our trap, the total number of atoms is iV = 2.2(0.3) x 10 5 . From 
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the measured trap frequencies, corrected for anharmonicity, we obtain u± = ^/uj x u y = 2ir x 
1696(10) Hz and u z = 2vrx72(5) Hz, so that a; = (u x u y u) z ) l/3 = 2?r x 592(14) Hz is the mean 
oscillation frequency. For these parameters, the typical Fermi temperature Tp = (3Ny/ 3 huj jks 
for a noninteracting gas is ~ 2.5 /xK, small compared to the final trap depth of Uo/kB = 35 fiK. 

Energy is precisely added to the trapped gas at fixed atom number by releasing the cloud 
from the trap and permitting it to expand for a short time t heat after which the gas is recaptured. 
As shown below, even for the strongly-interacting gas, the energy input is well-defined for 
very low initial temperatures, where both the equation of state and the expansion dynamics 
are known. During the times th eat used in the experiments, the axial size of the gas changes 
negligibly, while transverse dimensions expand by a factor b±(t heat ). Hence, the harmonic 
trapping potential energy in each of the two transverse directions increases by a factor b 2 ± (t heat ). 

The initial potential energy is readily determined at zero temperature. This follows from the 
equation of state of the gas, (1 + /3)e_F(x) + U tr a P {*) = where e_p(x) is the local 

Fermi energy, (3 is the unitary gas parameter (|7]|6]|3]|H]|3?3>, Utrap is the harmonic approximation 
to the trapping potential, and /z is the global chemical potential. The equation of state is 
equivalent to that of a harmonically trapped noninteracting gas of particles with an effective 
mass (0), which in our notation is M* = Mj (1 + (3), where M is the bare mass. Since the 
gas behaves as a harmonic oscillator, the mean potential energy is half of the total energy. As 
(3 < (8), M* > M, so that the effective oscillation frequencies and the chemical potential are 
simply scaled down, i.e., fi = k B T F y/l + (3 (|7]0). The total energy at zero temperature, which 
determines the energy scale, is therefore 

E = -Nfjio = -Nk B T F Jl^p. (1) 

For each direction, the initial potential energy at zero temperature is E /6. Then, the total 
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energy of the gas after heating is given by (31 ), 



E(t heat ) =f]E 



(2) 



Here, rj is a correction factor arising from the finite temperature of the gas prior to the energy 
input. For the noninteracting gas, r] nonint is determined at the lowest temperature T = 0.23 
from the energy for an ideal Fermi gas. For the strongly-interacting gas, where the initial 
temperature is very low and T = 0.04, we assume a Sommerfeld correction SJ2[ and obtain 
r] int ~ 1 + 27r 2 T 2 /3 ~ 1.01, which hardly affects the energy scale. 

The strongly-interacting gas exhibits hydrodynamic, anisotropic expansion ( 1 ), so that b± = 
b*[ is a hydrodynamic expansion factor (171135b. For the noninteracting gas, we use a ballistic 



expansion factor frf (i) = y 1 + (to±t) 2 . The temperature change during the expansion time 
theat < 460 [is must be very small, since the minimum value of T = 0.04 is measured by 
imaging the interacting cloud after 1 ms of expansion. Hence, the primary heating arises only 
after recapture and subsequent equilibration. 

Thermometry of strongly interacting Fermi gases is not well understood. By contrast, ther- 
mometry of noninteracting Fermi gases can be simply accomplished by fitting the spatial distri- 
bution of the cloud with a Thomas-Fermi (T-F) profile, which is a function of two parameters. 
We choose them to be the Fermi radius a x and the reduced temperature T/Tp. However, this 
method is only precise at temperatures well below 0.5 Tp, where a x and T/Tp are determined 
independently. At higher temperatures, where the Maxwell-Boltzmann limit is approached, 
such a fit determines only the product a 2 x T/Tp. We circumvent this problem by determining 
a x from a low temperature fit, and then hold it constant in the fits at all higher temperatures, 
enabling a one-parameter determination of the reduced temperature. 

For strongly interacting Fermi gases below the superfluid transition temperature T c , the spa- 
tial profile may contain normal and superfluid components (173b . However, experimentally and 
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theoretically, one finds that the spatial profile of a strongly interacting gas closely resembles 
a T-F distribution (I71l5?l). For this reason, T-F fits to the cloud profiles are commonly used 
to estimate the reduced temperature, which is often reported as T/T F , where T F is the Fermi 
temperature for a noninteracting gas. Analogous to the noninteracting case, we define an ex- 
perimental dimensionless temperature parameter T, which is to be determined by fitting the 
cloud profiles with a T-F distribution 051) . holding constant the Fermi radius of the interact- 
ing gas, a' x . Unlike two parameter fitting procedures, this single parameter method is stable. 
We find experimentally that T increases monotonically from the highly degenerate regime to 
the Maxwell-Boltzmann limit. This fitting procedure also leads us to define a natural reduced 
temperature scale 



which is consistent with our choice of fixed Fermi radius a' x , i.e., Muj x a x /2 = /i . At high 
temperatures, we must interpret T = T nat , to obtain the correct Maxwell-Boltzmann limit. At 
low temperatures, T ~ T nat yields an estimate of T/Tp. However, to determine the precise 
correspondence between T and the reduced temperature T/Tp which is input into theoretical 
models, one should perform the experimental fitting procedure with the theoretically generated 
density profiles as suggested and implemented by Chen et al., (l3oT) . 

The experimental fitting procedure measures T by first obtaining one dimensional, trans- 
verse spatial distributions n(x) from the column density by spatially integrating along the trap 
axial direction. Dividing by the total number of atoms per spin state, we obtain normalized spa- 
tial profiles. Then T is determined using the one parameter T-F fit method, yielding 0.04-2.15 
for the strongly-interacting gas and 0.2-1.1 for the noninteracting gas. 

The experimental energy scale Eq. |2] and the natural temperature scale Eq. |3] are deter- 
mined by calculating (3 from the measured Fermi radii for the interacting and noninteracting 



gas samples. The relation is given by a' x = a x (l + /5) 1//4 ©, where a x = ^2k B T F / '(Mu%) 
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is the radius for a noninteracting gas. To determine a' x , we measure the size of the cloud af- 
ter 1 ms of expansion, and scale it down by the known hydrodynamic expansion factor of 
b H (l ms) = 13.3 d7]Q3>. We then determine the Fermi radius a' x = 11.98 (iV/2) 1 / 6 /im/13.3 = 
0.901(0.021) (iV/2) 1 / 6 /i m - Using o x = 1.065 (iV/2) 1 / 6 /im for our trap parameters, yields 
[3 = —0.49(0.04) (f57b in reasonable agreement with the best current predictions, where (3 = 
-0.56 ©, and p = -0.545 §M- 

We now apply our energy input and thermometry methods to measure the heat capacity of 
an optically trapped Fermi gas, i.e., for different values of theat, we measure the temperature 
parameter T and calculate the total energy E(th ea t) / E from Eq. (2J To obtain high resolution 
data, 30-40 different heating times theat are chosen. The data for each of these heating times are 
acquired in a random order to minimize systematic error. Ten complete runs are taken through 
the entire random sequence. 

To test the method with a known system, we first measure the heat capacity for a noninteract- 
ing Fermi gas at 526 G. The gas is initially cooled to T = 0.23 (the lowest temperature we can 
achieve in this case) by 30 seconds of forced evaporation at 300 G as described previously dZJt. 
and then heated as described above. Fig. [T] shows the data (green dots) which represent the cal- 
culated E(t heat ) j E versus the measured value of T, for each t heat . For comparison, predictions 
for a noninteracting, trapped Fermi gas, Eid ea i(T) / Eid ea i(0) are shown as the red curve, where 
T = T/T F in this case. Here, the chemical potential and energy are calculated using a finite 
temperature Fermi distribution and the density of states for the trapped gas. We use the density 
of states for a gaussian potential well fi28k . rather than the harmonic oscillator approximation. 
This yields very good agreement at all temperatures. 

Next, we measure the heat capacity for the strongly interacting gas at 840 G. Here the 
gas is cooled to T = 0.04 and then heated. Fig. [l] shows the data (blue dots) which represent 
E (theat) I Eq versus the measured value of T, for each t hea t- Note that the temperature parameter 
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T varies by a factor of 50 and the total energy by a factor of 10. Remarkably, on a large scale 
plot, the data for the strongly interacting and noninteracting gases appear quite similar. 

A striking result is shown by plotting the data for the strongly interacting gas on a log — log 
scale. Fig. 12 shows that the data reveal a transition in behavior at T ~ 0.33, where the slope 
changes. Above T ~ 0.33, the data for strongly interacting data overlap closely with that of the 
noninteracting gas. Below T ~ 0.33, the data deviates significantly from noninteracting Fermi 
gas behavior. This transition may arise from changes in the behavior of the total energy and 
from changes in the spatial profile of the gas which serves as our thermometer. 

Insights into the microscopic structure of the strongly interacting gas can be obtained from 
the temperature scaling of the energy. Above the transition, for T > 0.33, we find that the data 
in Figures [Hand El are well fit by E{T) = ^TTp E ideal {T), with a constant (3 = -0.49. This 
suggests that T = T nat = T/(T F \/1 + (3) is a good approximation above the transition. Such 
scaling may be a manifestation of universal thermodynamics (Offl) . 

Below the transition, for T < 0.33, the gas may comprise several components, for example, 
a normal Fermi gas and both superfluid and noncondensed pairs, each contributing differently 
to the temperature scaling, as arises in a pseudogap model U6\\36h . For simplicitiy, we consider 
here a temperature scaling of the form E(T)/E = 1 + bT c . For sufficiently low temperature, 
one expects c = 2 for an ideal Fermi gas, c = 5/2 for a homogeneous noninteracting Bose gas, 
and c = 4 for a harmonically trapped Bose gas. The best fit (black line in Fig. 12) corresponds 
to c = 2.53(0.15) and b = 9.8(1.9). The x 2 P er degree of freedom for this fit is 1.4. We find 
that the parameters b and c are strongly correlated. Holding c = 5/2, we obtain b = 9.4(0.2). 
Fitting a quadratic temperature dependence yields b = 4.8(0.2), and a larger x 2 per degree of 
freedom of 5.2. A T 4 power law fit yields b = 63.8(4.1) and a x 2 P er degree of freedom of 
7.1. These results suggest that the gas is neither a normal Fermi gas nor a BEC of small weakly 
interacting molecules. 
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One can understand 5/2 power scaling at very low temperature as arising from thermal exci- 
tation of low energy bosons (fermion pairs) (18), where the fermionic contribution is exponen- 
tially suppressed by the superfluid gap. A simple picture of the 5/2 power scaling is that short 
wavelength thermal excitations increase the local kinetic energy of bound pairs without breaking 
them, yielding the density of states and energy for free particles in three dimensions. To make an 
estimate of b based on universal scaling, we assume that the bosons have a mass of 2M*, so that 
the density of states per unit volume for a locally homogeneous gas is (4M*/fi 2 ) 3 / 2 e 1 / 2 / (4-7T 2 ). 
The total energy is easily determined using a Bose distribution with zero chemical potential. 
Assuming that the pairing energy scale is large compared to k B T over most of the trap volume, 
the e integration is approximated from to oo. Multiplying the resulting energy density by 
the trap volume N/n, where n is the average density, we obtain E/E = 1 + bT^J, where 
b = (3/4)C(5/2)(2tt) 3 / 2 (315/512) = 9.75, close to the result 9.4(0.2) obtained from the fit. 

We estimate the transition temperature from the intersection point, T = 0.33(.02) (l37l) . of 
the power law fit and the scaled ideal gas prediction, Fig. |2j To extract a preliminary value of 
T c /Tp, we assume T = T nat near the transition temperature, and use the measured value of 
@ = -0.49. We then obtain T C /T F = y/TTfif = 0.24(.02) close to predictions for the 
superfluid transition temperature which have been made over the last decade (li #11391 17711?^) . 
The fractional change in the heat capacity C is estimated from the slope change in the fits to the 
data, assuming that the temperature calibration function is smooth near T c (36). In that case, 
(C> - C<)/C> = -0.48(0.03), where > (<) denotes above (below) T c . 

Recently, Q. Chen, J. Stajic and K. Levin have done a pseudogap model of a trapped, 
strongly interacting Fermi gas (36), and obtain both the energy and the spatial profile as a 
function of reduced temperature T/Tp, throughout the superfluid and normal region. The tem- 
perature scale T/Tp is calibrated to our T by fitting one dimensional T-F profiles to the the- 
oretical spatial distributions as described above, yielding a monotonic relation. The data of 
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Figures [T] and |2l are very well reproduced by the theory. 
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Figure 1: Total energy versus temperature. For each heating time th. ea u the temperature pa- 
rameter T is measured from the cloud profile, and the total energy E(theat) is calculated from 
Eq. El in units of the ground state energy E . Green circles: noninteracting Fermi gas data; 
Blue circles: strongly interacting Fermi gas data. Red curve: predicted energy versus reduced 
temperature for a noninteracting, trapped Fermi gas, E idea i(T) / E idea i(0). 
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Figure 2: Energy input versus temperature from Fig.[T]on a log — log scale. The strongly inter- 
acting Fermi gas shows a transition in behavior near T = 0.33. Green circles: noninteracting 
Fermi gas data; Blue circles: strongly interacting Fermi gas data. Red curve, prediction for a 
noninteracting, trapped Fermi gas. Black line, best fit power law 9.8 T 2 53 . Note the lowest 
temperature point (blue square) is not included in the fits, as it is constrained to lie on the red 
curve. 
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